Bound states in the continuum in a single-level Fano- Anderson model 
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Bound states in the continuum (BIC) are shown to exist in a single-level Fano- Anderson model with 
a colored interaction between the discrete state and a tight-binding continuum, which may describe 
mesoscopic electron or photon transport in a semi-infinite one-dimensional lattice. The existence of 
BIC is explained in the lattice realization as a boundary effect induced by lattice truncation. 
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I. INTRODUCTION. 



Since the pioneering proposal of the existence of iso- 
lated quantum mechanical bound states embedded in the 
continuum, made by von Neumann and Wigner [l| in 
the study of the one-particle Schrodinger equation with 
certain spatially oscillating attractive potentials, several 
theoretical and a few experimental studies have demon- 
strated the existence of "bound states in the contin- 
(BlCjin a wide range of different physical systems 
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221 ]. BIC may be found in certain atomic or 
molecular systems @, loUpl , such as hydrogen atom in a 
uniform magnetic field Jo], in semiconductor superlattice 
structures [H, H, H, [ll|, EH , in mesoscopic electron trans- 
port and quantum waveguides H, H , M , fiol. 12 ,[3, EE], and 
in quantum dot systems (l3l. ll7L Tl9L l20l. 21 1. The ques- 



uum 



tion of existence of BIC has been also addressed for the 
famous Fano- Anderson Hamiltonian [23|, [24], 25 1 (also re- 
ferred to as the Friedrichs-Lee Hamiltonian 26 . [27\, |28|). 
which model the process of quantum mechanical decay 
of an unstable localized state coupled with a contin- 
uum in different contexts such as atomic physics, quan- 
tum electrodynamics, solid state and high-energy physics 
(see, e.g., 0, [M H| HH HI). BIC in Fano- Anderson- 
like models are commonplace in case where several (i.e. 
more than one) discrete states are coupled to a com- 
mon continuum; a noteworthy example of this case in 
condensed-matter physics is provided, for instance, by 
quantum transport and scatte ring in dot molecules at- 
tached to leads [U El HI HI HI. In those sys- 
tems the existence of BIC is usually related to the de- 
struction of discrete-continuum decay channels via quan- 
tum interference through a typical trapping mechanism. 
The conditions for the existence of BIC for a general 
multi-level Fano-Anderson Hamiltonian have been re- 
cently stated in Ref . (33j ; in particular, a sufficient con- 
dition that ensures the non-existence of BIC has been 
demonstrated. In case of a single localized state em- 
bedded in and interacting with a continuum, the exis- 
tence of bound states (sometimes referred to as " dressed 
bound states") has been acknowledged on many occa- 
sions and related to threshold effects or to singularities 
or gaps in the density of states of the continuum (see, 
e.g., [H, [H H HE [H, H3, H, HI). However, such 



dressed bound states have usually an energy outside the 
continuum [25]: Effects such as fractional decay, popu- 
lation trapping and atom-photon bound states found in 
several physical models describing the decay of a single 
discrete level coupled to a continuum [3l|, H3, IH [Hj] 
are in fact most of the times related to the existence of 
dressed bound states with an energy outside the contin- 
uum. Conversely, BIC have been seldom encountered in 
the single-level Fano-Anderson model: Examples of BIC 
involving a single localized state coupled with a specially- 
structured continuum have been found in the study of 
certain exactly-solvable electrodynamic models of spon- 
taneous emission decay with a density of modes showing 
a point-like gap [H Hl[ , however physical realizations of 
such models were not proposed [3X1 ] . 
It is the aim of this work to present an exactly-solvable 
Fano-Anderson Hamiltonian with a single discrete level 
coupled to a continuum by a colored interaction which 
supports BIC. The proposed model is of physical rel- 
evance to certain condensed-matter and photonic sys- 
tems, since it may describe the charge transfer dynam- 
ics of adatoms to a semi-infinite one-dimensional lattice 
of quantum dots [42| or photon tunneling dynamics in 
semi-infinite optical waveguide arrays or coupled opti- 
cal resonators [43|, Sj, [45|, 5a S3, |48| within the tight- 
binding approximation. In these lattice realizations, the 
tight-binding analysis allows one to simply explain the 
appearance of BIC as a boundary effect due to trunca- 
tion of the lattice. 

The paper is organized as follows. In Sec. II, the Fano- 
Anderson model describing the coupling of a single lo- 
calized state with a continuum is briefly reviewed, and 
the conditions for the existence of bound states either 
inside or outside the continuum are presented. Section 
III proposes an exactly-solvable model with colored in- 
teraction which supports BIC; exact analytical results for 
the fractional decay induced by BIC are also presented. 
Finally, in Sec. IV a lattice realization of the model is 
presented, together with a simple physical explanation 
for the existence of BIC. 
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II. BOUND STATES AND DECAY DYNAMICS 
FOR THE SINGLE-LEVEL FANO-ANDERSON 
MODEL 

A. Basic Model 

The starting point of the analysis is provided by a stan- 
dard Fano- Anderson model describing the interaction of 
a discrete state |a), of energy huj a , with a continuum 
described by a set of continuous states |fc) with energy 
Huj(k) (see, for instance, [H, HSIllI)- The Hamiltonian 
of the interacting system can be written as H = Hq + V, 
where 

H = hw a \a){a\ + J dk Mj(k)\k)(k\ (1) 

is the Hamiltonian of the non-interacting discrete and 
continuous states, and 

V = hJ dk[v{k)\a){k\+v*{k)\k){a\] (2) 

is the interaction part. Normalization has been assumed 
such that (a\a) = 1, (a\k) = and (k'\k) = S(k - k'). If 
we expand the wave function \ip) of the system as \tp) = 
c a {t)\a) + J dkc(k,t)\k), the expansion coefficients c a (t) 
and c{k 1 t) satisfy the coupled-mode equations 

ic a (t) = uj a c a + J dkv(k)c(k,t) (3) 
ic(k,t) = uj{k)c{k,t) + v*(k)c a (t), (4) 

where the dot indicates the derivative with respect to 
time. Typically, we assume that the frequency uj(k) 
of continuous states spans a hnite interval (a band) 
lui < u> < 0J2 for the allowed values of the continuous 
variable k, and that the frequency ui a of the discrete level 
is embedded in the continuum, i.e. that lo\ < uj a < lu2- 



B. Bound states 

The eigenstates \"4>e) of H corresponding to the eigen- 
value E — Ml are obtained from the eigenvalue equation 
H\4>e) = Kl\%j}E)- After introduction of the density of 
states p(cu) — dk/dui, from Eqs.(l) and (2) it follows that 
the eigenfrequencies il are found as the eigenvalues of the 
system 

VLc a = Ul a C a + / doj-y/ p(oj)v(u))c(w) (5) 
J ui 

Qc(u) — ujc(uj) + v p(uj)v* (uj)c a , (6) 

where we have set c(u) = \J p(uj)c(ui). As the continu- 
ous spectrum of H is the same as that of Hq, isolated 
eigenvalues corresponding to bound states may or may 
not occur for H. The energy E corresponds to a bound 



state of H provided that \iPe) is square integrable. This 
implies 

li^|| 2 = M 2 + / ^|5H| 2 <oo. (7) 

We have to distinguish two cases. 

(i) Bound states outside the continuum. This is the 
most common case, which has been studied on several 
occasions. A bound state with energy Ml outside the 
continuum exists provided that a root of the equation 

n - uj a = A(n) (8) 

can be found outside the band (wi,^), where we have 
set 

A(0)= r^(»)M">l a . (9) 
•V - w 

The conditions for the existence of bound states of 
this kind have been extensively investigated in the 
literature (see, e.g., [U, [H, For instance, Eq.(8) 

admits always two solutions outside the interval (oji,^) 
whenever p(ui)\v(ui)\ 2 does not vanish at the edge of 
the band, since in this case A(f2) diverges to =Foo as 

— > UJj 2 - 

(ii) Bound states inside the continuum. As shown 
in Kef. ,'!.''> . a bound state at frequency $1 inside the 
continuum may exist provided that the following two 
conditions are simultaneously satisfied 

\v(n)\ 2 P (n) = o , n-u a = A(Q) (10) 

Additionally, v(u)y/p(uj) should vanish as u> — > fl at least 
as ~ (lo—D,) in order to ensure a finite norm [Eq.(7)]. The 
first equation in (10) can be satisfied for either p(O) = 
or v(Q) = 0. The former case, which corresponds to a 
point-like gap in the density of states inside the band, 
has been previously considered for some special density 
of state profiles [4CJ, [H| , which however do not seem to 
have simple physical realizations [3l|. The latter case, 
u(O) = 0, implies that the discrete state \a) does not 
interact with the continuous state of frequency f2, and 
thus implies a "colored" interaction profile v(u>) with one 
zero at lo = f2. However, at such a frequency the addi- 
tional condition uj a — Q — A(fi) must be simultaneously 
satisfied, which means that BIC may exist solely at a 
prescribed energy fku a of the level \a). 

C. Decay dynamics 

Consider now the decay dynamics of the unstable 
state | a) embedded in the continuum. This corre- 
sponds to solve Eqs.(3) and (4) with the initial condi- 
tions c a (0) = 1 and c(fc,0) = 0, which can be done 
by e.g. a Laplace transformation or a Green's function 
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analysis (see, for instance, [25|, [29|, [3l|)- Indicating by 
c a (s) = / °° dtc a (t) exp(— st) [Re(s) > 0] the Laplace 
transform of c a (t), from Eqs.(3) and (4) one obtains 



discrete state with a continuum, which admits of BIC. 
Precisely, we assume the following colored interaction 
function 



c a {s) 



IS - 0J a 



and then, after inversion 



Ca(t) 



+ +ioo 



ds- 



exp(si) 



2n J +_ ioo ~ is - w a - T,(s) 
where S(s) is the self-energy, given by 

l«(*)| 2 



E( S ) 



dfc- 



w(fc) 



(ii) 



(12) 



(13) 



Possible poles on the imaginary axis of c a (s) correspond 
to bound states of H and are responsible for fractional 
decay of the amplitude c a (t). In fact, using the property 



S(s = -iu ± 0+) = A(lo) =f i7rp(o;)|w(w)| 2 



with 



J LOl 



OJ — UJ 



(14) 



(15) 



the poles s p = —iQ of c Q (s) satisfy the conditions Q — 
u> a = A(O) and p(f2)|?;(f2)| 2 = 0, i.e. they are located 
in correspondence of the bound states (either outside or 
inside the continuum) of H. In absence of poles (i.e. 
of bound states of H), c a (t) decays to zero, whereas in 
presence of poles Eq.(12) can be written as the sum of a 
contour (decaying) integral plus the (non-decaying) pole 
contributions. 



v{k) 



K a sm(nok) 



(16) 



and the following tight-binding dispersion curve for the 
band of continuous states 



i{k) = — 2kq cos k 



(17) 



where < k < n, no and n a are positive real parameters, 
and no is a positive and nonvanishing integer number, 
i.e. no = 1,2,3, .... A physical realization of this model 
will be described in the next section. We note that this 
model is a generalization of the well-known tight-binding 
Fano- Anderson model with a constant interaction cou- 
pling v(k) = const, which is known to show bound states 
outside the continuum [25| . Note also that for this band 
model the density of states, given by 



1 



dk 

P^' = IT = ~ rr-> ^ 



(18) 



shows van-Hove singularities at the band edges u) — 
±2ko- For a constant coupling 25], these singularities 
are responsible for the existence of two bound states out- 
side the band from either sides. However, for the colored 
coupling considered in our case [Eq.(16)] one has 



III. BIC IN A FANO-ANDERSON MODEL 
WITH A COLORED INTERACTION: AN 
EXACTLY-SOLVABLE MODEL 

In this section we present an exactly-solvable Fano- 
Anderson model, describing the interaction of a single 



G(u) = p(w)|«H| s 



2k 



npacos 



(19) 

which vanishes at the band edge. From Eqs.(13), (16) 
and (17), the self-energy E(s) can be calculated in an 
exact form and reads 



S( S ) 



2k| 

7T 



dk- 



sin 2 (nofc) 



IK' 



is - 



2kq cos k 



+ a4 



i - 



4«g 
2k 



2n 



(20) 



Using Eq.(14), the following expression for A(u>) = 
Re[E(s = —iuj ± + )] can be then derived 
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A(w) = < 



sin 



1 



2k 



2n 



^/co 2 -4reg 



1 



2no acos (^) 



2k () 



2n 



, a; < — 2ko 

, — 2k < w < 2k 



(21) 
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FIG. 1: (color online) Behavior of G(lo) = p(ui)\v(uj)\ 2 (dot- 
ted curves) and A(w) (solid curves), normalized to kq, versus 
normalized frequency cj/kq for (K a /Ko) = 0.2 and for increas- 
ing values of integer no. 



The behavior of G(u) = P (uj)\v(uj)\ 2 [Eq.(19)] and 
A(w) [Eq.(21)] for increasing values of no is shown in 
Fig.l. Note the oscillatory behavior of both G(u) and 
A(w), with the existence of (2no — 1) zeros of A(u) 
at loi = — 2«o cos[/7r/(2n )] (I = l,2,...,2n — 1) and 
of (no + 1) zeros of G(u>) at uj m — — 2kq cos(mir/no) 
(to = 0, 1, 2, no). We can then specialize the general 
results of Sec. II to the present model. 

(i) Bound states outside the continuum. At most 
two bound states at frequency Q outside the band 
(— 2/«o,2ko) from either sides may exist. Precisely, a 
bound state at frequency f2 > 2no exists provided that 
2kq — uj a < A(2«o), i-e. ui a > 2kq — n^no/no, whereas 
a bound state at frequency £1 < —2k does exist for 
— 2kq — uj a > A(— 2kq), i.e. for ui a < —2kq + k^toj/ko- 
Therefore, if the frequency uj a of the discrete level lies 
inside the interval 



1 + #<^<1 
2kq 2ko 



2kq 



(22) 



no bound states outside the continuum do exist. For a 
given value of no, Eq.(22) is satisfied for a sufficiently 



small value of the coupling K a j kq and provided that the 
frequency u a is not too close to the band edges. In the 
following, it will be assumed that no bound states exist 
outside the continuum. 

(ii) Bound states inside the continuum. Following 
the results of Sec. II. B, one BIC for the model expressed 
by Eqs.(16) and (17) does exist at = ui a for any 
n o > 2, provided that the frequency ui a of the discrete 
level assumes one of the following (no — 1) allowed 
values: 

uj a — — 2kq cos(m7r/no) (to = 1, 2, no — 1). (23) 



(iii) Decay dynamics and fractional decay due to BIC. 
Suppose that H admits of one BIC but no bound states 
outside the continuum, i.e. that Eqs.(22) and (23) are 
simultaneously satisfied. The decay law for c Q (t) is given 
by the inverse Laplace transform Eq.(12). The Bromwich 
integration path in Eq.(12) can be deformed into the con- 
tour a shown in Fig. 2, where c a (s) is always calculated on 
the first Riemannian sheet. The integral then comprises 
the pole contribution at s p = —iuj ai which arises from the 
semi-circles surrounding the pole, and the principal- value 
integral of [1 / (2ni)]c a (s — — iuj±0 + ) exp(-iujt) along the 
interval — 2ko < uj < 2kq of the imaginary axis from the 
two sides Re(s) = ±0 + of the branch cut, i.e. 



Ca(^) — Cpole(t) ~t~ CcZecay(^)- 



(24) 



Using Eqs.(12), (14) (19) and (21), after some straight- 
forward calculations one then obtains: 



exp(-iuj a t) 



n ( Kg \ 

2 \«oj 









t 


i - 






\2k-o J 



-1/2 



(25) 



for the pole contribution (non-decaying term), and 
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1 ' ^ 2 



2ir \ Kq 



Cdecayit) = — ( ^ ) V \ dk X (26) 



sin 2 (no£;) exp(2iKoi cos k) 



[u> a / (2k ) + cos k — (K a /2n ) 2 sin(2n fc)/ sinfc] 2 + [k 2 /(2kq)] 2 sin 4 (nofc)/ sin 2 A; 



Re(s) t 



-2iK„ 



9- 



2nc Im M 



FIG. 2: Deformation of the Bromwich path B (dotted line) 
for inverse Laplace transformation. The bold solid curve is the 
branch cut, whereas the deformed path is represented by the 
solid closed contour a surrounding the branch cut. The BIC 
corresponds to the pole s p = — iui a of c a (s) on the imaginary 
axis internal to the branch cut, which is surrounded by two 
semi-circles whose radius R tends to zero. 



for the decay term. The existence of a BIC is thus re- 
sponsible for a fractional decay of the amplitude c a (t). 
Such a fractional decay is different from the most com- 
mon one encountered in other single-level Fano- Anderson 
models (see, e.g., [3!|) since in our model the fractional 
decay is due to the existence of a BIC. Note also that, 
if no bound states exist, i.e. if Eq.(22) is satisfied but 
one of the resonance conditions (23) is not satisfied, the 
amplitude c a (t) fully decays toward zero. In this case, 
one simply has c a (i) = c deC ay(t), where c de cay(t) is given 
again by Eq.(26) in which the principal value of the in- 
tegral may be omitted. It is worth considering the limit 
tiq — ► oo for the model expressed by Eqs.(16) and (17). 
In this case, in the integral on the second right term of 
Eq.(20) the rapidly-oscillating function sin 2 (nofc) can be 
replaced by its cycle-averaged value 1/2, i.e. one can 
approximately write 



k 2 r 1 

£(s) ~ ^ / dk- 



-IK„ 



/0 is + 2k cos k ^/s 2 + 4k 2 , 
which from Eq.(14) yields 

G(u) 



and 



uj < —2kq 



(27) 



(28) 



A(w) = { -2k < lu < 2k (29) 

Note that the above expressions for G(u>) and A(w) cor- 
respond to the limit of a tight-binding Fano- Anderson 



model with an uncolored interaction (see, e.g., Ref . [251] . 
pp. 283-285), i.e. to a flat interaction function v(k) ~ 
Ka/y/Tr. In this case it is knwon [25[ that BIC do not exist 
at any value of tu a , whereas two bound states outside the 
continuum are always found. The physical explanation 
of the disappearance of BIC in the uq — > oo limit will be 
discussed in the next section. 



IV. A TIGHT-BINDING LATTICE 
REALIZATION OF THE COLORED 
FANO-ANDERSON MODEL 

In this section we propose a simple and noteworthy 
physical realization of the Fano- Anderson model with col- 
ored interaction discussed in the previous section, which 
may describe either electron or photon transport phe- 
nomena in condensed-matter or photonic tight-binding 
lattices. Since k varies in the range < k < 7r, we can 
expand c(fc,t) as a Fourier series of sine terms solely ac- 
cording to 



c(k,t) 



(30) 



c n (t) sin(nfc) 

n=l 

where the time-dependent coefficients c n are given by 

c n (t) = —\ — [ dk c(fc, t) sin(nfc). (31) 
V 7T Jo 

Taking into account that 

/ dk sin(nfc) sin(mfc) = —5 n ,m (32) 
Jo 2 

(n,m > 1), from Eqs.(4), (16), (17) (30) and (31) the 
equations of motion for the coefficients c n can be readily 
derived and read 

ic n = -K (c n+ i + Cn-l) - K a C a 5n,n (n > 2) (33) 

ici = -K Q c 2 - K a c a S n „ t i. (34) 
The equation for c a [Eq.(3)] then reads 

In their present form, Eqs. (33-35) can be derived from 
the tight-binding Hamiltonian 



H T b = -h2_^K Q (\n)(n + l\ + \n + l)(n\) + huj a \a)(a\ + 

n=l 

-hK a (\a)(n \ + \n )(a\) (36) 
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FIG. 3: (color online) Schematic of a localized state \a) cou- 
pled to a tight-binding semi-infinite lattice which realizes the 
Hamiltonian model (36). ko is the hopping amplitude be- 
tween adjacent sites in the lattice, whereas K a is the hopping 
amplitude between the localized state \a) and the site \no) of 
the lattice. 



which describes the interaction of the localized state \a) 
with the no-th site of a semi-infinite one-dimensional 
tight-binding lattice in the nearest-neighbor approxima- 
tion (see Fig. 3). The tight-binding model expressed by 
Eq.(36) has been used to study transport phenomena in 
different physical systems, including photon tunneling 
dynamics in coupled optical waveguides [isj or in cou- 
pled photonic cavities [43|, [44], HE [47j , charge transfer of 
adatoms to a one-dimensional lattice of quantum dots 
[42| , or decay of the polarization in spin chains [4!| . For 
instance, Eq.(36) may describe charge transfer between 
an adatom localized state and a one-dimensional mini- 
band associated with a quantum dot array, the adatom 
being attached to the semiconductor quantum-dot array 
surface [12] • It should be noted that these previous mod- 
els considered cither a semi-infinite tight-binding chain 
with a boundary defect [H, [49[ corresponding to the spe- 
cial case no = 1. where however no BIC exist, or to an 
infinite array [25ll42T|. i.e. to uq — > 00 corresponding to a 
non-colored interaction, where again no BIC exist. How- 
ever, as shown in the previous section, for a finite value 
of no larger than one, i.e. by considering a semi-infinite 
lattice in which a defect state interacts with a lattice site 
\rio) near (but not at) the boundary, BIC at certain fre- 
quencies w a do exist according to Eq.(23). In the tight- 
binding representation (36) the existence of BIC has a 
simple physical explanation which is related to a bound- 
ary effect of the semi-infinite lattice: BIC correspond to 
localized states in the lattice with c n = for n > no- In 
fact, let us look for a solution to Eqs. (33-35) of the form 
c„ = c n exp(— iflt), c a = c a exp(— i£lt), with c n = for 
n > no- From Eq.(35) and Eq.(33) with n = no it then 
follows that Q = u) a and c a = — (ko/k q )c„ _i, whereas 
after setting c = {c\, C2, c„ _i) T from Eq.(34) and 
from Eq.(33) with n < ng — 1 one obtains that c and £1 
are the eigenvectors and corresponding eigenvalues of the 
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Site Number n 

FIG. 4: (a) Decay dynamics of the amplitude |c a (f)| as ob- 
tained by numerical analysis of Eqs. (33-35) for parameter val- 
ues Ko = 1, K a = 0.2, no = 12, and for u> a = 0.15, correspond- 
ing to the absence of BIC. In (b) the temporal evolution of 
the amplitudes |c n (f)| at the lattice sites is also shown on a 
grey-scale plot. 
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(37) 

i.e. Aic = uo a c Diagonalization of the matrix M. yields 
for the eigenvalues the following expression 

Q m — — 2kq cos(m,7r/no) (rn = 1, 2, no — 1) (38) 

with corresponding eigenvectors = sin(m7rn/no). 

Note that the values of Vl m given by Eq.(38) are precisely 
the resonance frequencies for the existence of BIC found 
in the previous section [see Eq.(23)]. Therefore in the 
tight-binding realization of the colored Fano- Anderson 
model BIC arise due to a trapping effect which localizes 
the excitation between the boundary of the semi-infinite 
chain and the |no — l)-th site of the chain. The cou- 
pling of the state \a) with the lattice site |no) allows 
for the vanishing of the amplitudes c„ at lattice sites 
n > no via quantum destructive interference. It should 
be noted that similar trapping mechanisms supporting 
BIC in tight-binding models have been recently found in 
trip le or quadruple dot molecules connected to two leads 
[2Q , [22j], which are modeled as two semi- infinite tight- 
binding lattices. In these models, BIC correspond to van- 
ishing of the wave function at the sites in the molecule 
in contact with the leads, i.e. BIC are fully localized in 
the molecule sites but not in the leads. Conversely, the 
present tight-binding model [Eq.(36)] involves solely one 
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FIG. 5: Same as Fig. 4, but for uj a — 0, corresponding to the 
existence of one BIC. 

localized state side-coupled to a semi-infinite lattice, and 
therefore BIC can not simply correspond to a decoupling 
of the localized state with the lattice. This is clearly 
demonstrated by the fact that the wave function of a 
BIC for the model expressed by Eq.(36) is non- vanishing 
even in a portion of the lattice (from the boundary site 
|1) to the site \n — 1)). Additionally, BIC cease to exist 
as riQ — > oo, i.e. lattice truncation is essential to sustain 
BIC. 

We checked the existence of BIC induced by this trapping 
mechanism by a direct numerical analysis of the coupled 
mode equations (33-35) using a fourth-order variable-step 
Runge-Kutta method with smoothly absorbing bound- 
ary conditions at the right boundary of the lattice to 
avoid spurious reflections due to truncation of Eq.(33). 
As an example, Figs. 4 and 5 show the decay dynamics 
of |c Q (i)| for parameter values corresponding to a com- 
plete decay (Fig. 4), i.e. to the absence of a BIC, and to 
a fractional decay (Fig. 5) related to the existence of a 
BIC. In the figures, the dynamical evolution of the lat- 



tice site amplitudes |c„(£)| is also depicted on a gray- 
scale plot, showing either a diffusion [Fig. 4(b)] or a local- 
ization [Fig. 5(b)] of the excitation transferred from the 
localized state \a) to the lattice site |no). We checked 
that the numerically-computed decay law for c a (t) ex- 
actly reproduces the curve predicted by the analytical 
decay law Eqs. (24-26). Note that, as the resonance con- 
dition (26) for the existence of BIC is satisfied, fractional 
decay of c a is attained [Fig. 5(a)], which clearly corre- 
sponds to trapping of the excitation at the lattice sites 
|1), 1 2), |no — 1) with a destructive interference of site 
excitation for n > no [Fig. 5(b)]. Conversely, for a fre- 
quency uj a which does not satisfy the resonance condi- 
tion (23) for some integer m, the amplitude c a (t) fully 
decays toward zero [Fig. 4(a)] and the excitation trans- 
ferred to the lattice diffuses along the lattice without be- 
ing trapped [Fig. 4(b)]. 

V. CONCLUSIONS 

In this work an exactly-solvable single-level Fano- 
Anderson model which admits of bound states inside 
the continuum has been proposed, and its relevance to 
tight-binding lattice models generally adopted to study 
electron or photon transport phenomena in condensed- 
matter or photonic systems has been discussed. As pre- 
viously proposed models supporting BIC in single-level 
Fano- Anderson models require point-like gaps in the den- 
sity of states [3l|, ES EJ- a condition which does not seem 
to have simple physical realizations [HJ- in the present 
work it has been shown that BIC can exist within a tight- 
binding continuum without point-like gaps provided that 
the interaction of the localized state with the continuum 
shows a colored interaction with zero-points. A lattice re- 
alization for such a colored Fano- Anderson model, which 
may be of relevance to model photon or electron trans- 
port in certain photonic or condensed-matter systems 
[H El EI HE EE SI HI HI , has been proposed, and a 
simple physical explanation of the existence of BIC as a 
trapping effect sustained by lattice truncation has been 
highlighted. 
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